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Abstract 

We extend a previous bispectrum analysis of the COBE-DMR 4 
year maps, allowing for the presence of correlations between all possible 
angular scales. We find that the non-Gaussian signal found earlier for 
bispectrum components including adjacent modes does not extend to 
triplets of modes with larger separations. Indeed for all separations 
Al > 1 we find that the COBE-DMR data is very Gaussian. The 
implication seems to be that the previously detected non-Gaussian 
scale-scale correlation falls off very quickly with mode separation. 



1 Introduction 

The issue of whether or not the primordial fluctuations are Gaussian is 
of fundamental importance in theories of structure formation. Gaussian 
large-scale primordial fluctuations are believed to be the hallmark of the 
simplest inflationary scenarios. In these theories vacuum quantum fluc- 
tuation, with Gaussian statistics, are stretched by a period of exponen- 
tial expansion to give rise to the large scale fluctuations we observe to- 
day. In contrast, topological defect theories, which are also possible candi- 
dates for structure formation |Contaldi 00, Pontaldi et al 99a ], are believed 
to produce non-Gaussian density fields. Non-Gaussianity is also believed 
to be associated with some non-minimal inflationary models [Salopek 1992, 
Linde and Mukhanov 1997 , Contaldi et al 99~b| , such as those generating isocur- 



vature fluctuations [ Peebles 97 ]. 

Several recent [ Ferreira, Magueijo &; Gorski (1998) , [Pando et al (1998} 



Magueijo (2000) ] papers have shown indications of non-Gaussianity in the 
cosmic microwave background (CMB) temperature fluctuations, as mea- 
sured by the COBE-DMR experiment. The first of these detections 
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[ Ferreira, Magueijo & Gorski (1998) ] was later found to be due to an experi- 
mental systematic [ Banday et al 9S |, and the second [ Pando et al (1998) ] to 
an error of method |Barrciro et al OC]. The third claim [ Magueijo (2000)' ], 
however, remains unassailable. In that work a bispectrum analysis was 
carried out, using adjacent scales {£ — 1,1, 1 + 1} as components of the bis- 
pectrum. A significant deviation from Gaussianity was found, with a high 
confidence level. 

An obvious question is whether this non-Gaussianity extends to compo- 
nents correlating scales with a larger separation, say {£ — A£, £,£ + A£}, with 
A£ > 1. The purpose of this paper is to answer this question. As we shall 
see the non-Gaussian signal found in | Magueijo (2000) | does not appear to 
extend to larger separations A£. 



2 Method 

In this letter we examine the possibility of deviations from Gaussianity in 
terms of the bispectrum, which should be zero for a Gaussian process. Let 
us expand the temperature fluctuations -^-(n) in Spherical Harmonic func- 
tions: 

~jr( n ) =^2 a £mY em (n) (1) 

tm 

A Gaussian probability distribution function is completely described by the 
first two moments, (a^ m ) and (a£ m agi m i) i.e. all higher moments can be 
obtained from the mean value and the variance. This does not hold for Non- 
Gaussian functions. The Bispectrum, a rotationally invariant cubic form 
associated with the third moment, (a.i m &i' m' m") •, is therefore a possible 
test of Gaussianity. It is given by 
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where the (...) is the Wigner 3 J symbol. In |Magueijo (2000) ] correlations 



between multipoles separated by A£ = 1, {£ — l£i + 1} were investigated. 
In this letter we complete this work by using the other possible separa- 
tions, Al = 2,3,4,5, and 6, which fit within the signal dominated range of 
I probed by the COBE-DMR experiment. Selection rules require ^1+^2 + ^3 



2 



to be even, |£$ — £j\ < £f. < i{ + £j and mi + m,2 + m.3 = 0. From these 
constraints we see that a suitable chain of correlators is Af~ £ = Bi-Mit+Aii 
with £ even. 



We then follow precisely the same procedure documented in | Magueijo (2000) | 



and [Ferreira, Magueijo & Gorski (1998)], apart from considering the more 



general ratio 

Jt l = ~. ; " , , ~ (3) 

{Ct-MYHCeYHCt+M) 1 / 2 

where C| = ^pr Em l a ^n| 2 - This quantity is dimensionless, and is in- 
variant under rotations and parity. The same data is examined; the in- 
verse noise- variance- weighted average maps of the 53A,53B,90A and 90B 
COBE-DMK channels, with monopole and dipole removed, at resolution 6, 
in galactic and ecliptic pixelization. We use the extended galactic cut of 
[|Banday et al 1997 1, and [Bennet et al 1996] to remove most of the emission 



from the plane of the Galaxy. To estimate the <//^s we set the value of the 
pixels within the galactic cut to and the monopole and dipole of the cut 
map to zero. We then integrate the map multiplied with spherical harmonics 
to obtain the estimates of the ae m s and apply equations |2| and ||. 

The observed J^ e s are to be compared with their distributions P(J^ e ) as 
inferred from Monte Carlo simulations in which Gaussian maps are subject 
to DMR noise and galactic cut. In simulating DMR noise we take into ac- 



count the full noise covariance matrix, as described in | Lineweaver et al 1994 ] 



This includes correlations between pixels 60° degrees apart. 



3 Results and conclusions 

The comparison between P(J^) and the observed for various A£ is dis- 
played in figs, [l] - |||. The observed values agree well with the probability 
distributions from the Monte Carlo distributions: implying that as far as 
A£ > 1 bispectrum components are concerned the data supports Gaussian- 
ity. 

To quantify just how well the observations agree with Gaussianity, we 
do a goodness of fit statistic on the data. We calculate the x 2 value for 
the different £ separations, and compare with the respective distributions 
F{\ 2 ). The rationale behind considering separate x 2 was explained in 
[fVlagueijo (2000)|1 P1, and for the distributions P{f l ) found (which are very 



1 This is to be contrasted with more dubious philosophy towar ds combining separate 



tests expressed in [Bromley and Tegmark 99, Phillips & Kogut 0C] 
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close to Gaussians) it makes sense to use the standard x 2 '- 

1 N ( \ 
X 2 = ^E 1 ^ ( 4 ) 

i=l 1 

where N is the number of different observables and /ij and af are the ex- 
pectation value and the variance of the ith observable respectively. This 



reduces to the expression proposed in [ Ferreira, Magueijo fc Gorski (1998) ] 
in our case. The F(x 2 ) are then calculated through further Monte Carlo 
simulations, so that possible correlations between the different I may be 
taken into account. The x 2 distributions and corresponding observed values 
are shown in fig ||. 

As we can see x 2 ~ 1 f° r an separations Al > 1, with deviations well 
within the spread expected from F(x 2 )- This is to be contrasted with the 
results obtained for Al = 1, that is bispectrum components measuring cor- 
relations between adjacent modes. For these x 2 ^ 1> meaning that the 
observed Je do not display the scatter around zero expected from a Gaus- 
sian distribution. We now found that this distinctive non-Gaussian signal 
does not extend to bispectrum components with Al > 1. The implica- 
tion seems to be that the inter-mode correlations suggested by the findings 



of [ Magueijo (2000) | fall off quickly with the scale separation, something 
which conforms to theoretical prejudice. 

In a recent paper [ Phillips Kogut 00) ] a variety of components of the 



bispectrum was studied in the context of a semi-analytical texture model. 
Such a model is at the very least oversimplified, and probably has little to 
do with real-life textures. In any case that work stresses the importance of 
analysing off-diagonal components, as we have done in this paper. 
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Figure 1: The vertical thick dashed lines represent the values of the observed 
J^ =2 , for different values of £. The solid lines are the corresponding proba- 
bility distribution functions of the jf^ =4 s for a Gaussian sky with extended 
galactic cut and DMR noise, as inferred from 10000 realizations. 
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Figure 2: The vertical thick dashed lines represent the values of the observed 
J^ =3 , for different values of £. The solid lines are the corresponding proba- 
bility distribution functions of the jf^ =4 s for a Gaussian sky with extended 
galactic cut and DMR noise, as inferred from 10000 realizations. 
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Figure 3: The vertical thick dashed lines represent the values of the observed 
J^ =4 , for different values of £. The solid lines are the corresponding proba- 
bility distribution functions of the jf^ =4 s for a Gaussian sky with extended 
galactic cut and DMR noise, as inferred from 10000 realizations. 
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Figure 4: The vertical thick dashed lines represent the values of the ob- 
served J^s, A£ = 5, and 6 for different values of i. The solid lines are 
the corresponding probability distribution functions of the J^ e s for a Gaus- 
sian sky with extended galactic cut and DMR noise, as inferred from 10000 
realizations. 
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Figure 5: The dashed lines are observed values of x 2 from the COBE data for 
different I separations A£ = 2, 3, 4 and 5. The solid lines are corresponding 
probability distribution functions for % 2 for a Gaussian sky with extended 
galactic cut and DMR noise, as inferred from 10000 realisations. 
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